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The feasibility of high-fidelity single-qubit operations of a hole spin in a quantum dot molecule by 
electric g tensor control is demonstrated. Apart from a constant external magnetic field the proposed 
scheme allows for an exclusively electric control of the hole spin. Realistic electric gate bias profiles 
are identified for various qubit operations using process-tomography-based optimal control. They 
are shown to be remarkably robust against decoherence and dissipation arising from the interaction 
of the hole with host-lattice nuclear spins and phonons, with a fidelity loss of ~ 1 percent for gate 
operation times of w 10 ns. Spin-echo experiments for the hole spin are modeled to explore dephasing 
mechanisms and the role of pulse-timing imperfections on the gate fidelity loss is discussed. 



I. INTRODUCTION 

We propose and study feasibility of all-electric con- 
trol of a qubit realization based on the hole spin in a 
quantum dot molecule. This system offers two main ad- 
vantages over electron-spin-based realizations: the use of 
hole spins increases the dephasing time associated with 
the interaction with nuclear spins by about an order of 
magnitude and allows for an efficient g tensor control, 
thereby facilitating essentially all-electric control of the 
qubit. 

Spin-based quantum bit realizations in semiconductor 
quantum dots have gained wide interest in the quantum 
computing community due to their potential regarding 
scalability and their long relaxation times. Most work, 
both theoretical and experimental, has focused on elec- 
tron spin qubits, which have spin relaxation times rang- 
ing from several milliseconds 1 up to secondsP However, 
it has turned out that this type of qubit is prone to 
decoherence processes due to the interaction with sur- 
rounding nuclear spins of the host lattice!^ The Fermi 
contact hyperfine interaction ultimately limits the inho- 
mogeneous dephasing time to T 2 * e ss 10 nsP Several tech- 
niques to circumvent this problem have been devised, e.g. 
nuclear state preparationpnS fabrication of silicon-based 
semiconductor heterostructures with zero nuclear mag- 
netic mo ment p or employing hole spins instead of elec- 
tron spins EEED The p-type symmetry of the hole Bloch- 
function leads to a cancelation of the Fermi contact hy- 
perfine interaction. However, it has been shown that the 
dipole-dipole hyperfine interaction and the coupling of 
the hole orbital angular momentum to the nuclear spins 
cannot be neglectedP^HH] These lead to an inhomoge- 
neous dephasing time T£ h one order of magnitude longer 
than that of electrons, which is in good agreement with 
recent experiments.^ 

Basic spin qubit implementation schemes usually re- 
quire local magnetic fields for individual spin manipula- 



tions.^ Another technique, which has been proposed and 
tested recently, is to exploit the tunability of electron- 
spin and hole-spin g tensors by electric fields rather than 
employing control by local magnetic fieldsP^Hin] if an ex _ 

ternal magnetic field B is present, the spin experiences 
an effective field g(E) ■ B that can be changed locally 
by means of the electric field E. The dependence of the 
g tensor on E is particularly pronounced for hole spins 
in vertically stacked quantum dot molecules since the lo- 
calization of the hole wave-function is very sensitive to 
externally applied electric fieldsP^The goal of this paper 
is to show that such a system provides a fast and efficient 
universal gate for simple pulse shapes of 10 ns duration. 

In Sec.|Hj we describe the computational basis states of 
the qubit, the corresponding Hamiltonian and the mech- 
anisms that lead to non-unitary dynamics. The equation 
of motion for the qubit, including decoherence and relax- 
ation dynamics, is given. In Sec. Ill we apply optimal 



control theory in order to find simple electric pulse shapes 
that execute the Hadamard gate, as well as a ir/2 and a 
7r pulse. We find that these transformations can be im- 
plemented with remarkably low fidelity losses of ~ 1%. 
In addition, we show how spin-echo experiments can be 
performed on the hole spin to explore decoherence in the 
system. We also discuss quantitatively the effect of pulse- 
timing imperfections. We conclude with a summary of 
the present work. The model for the interaction of the 
nuclear spins with the hole is detailed in Appendix |A"| 
and [0 



II. THEORY 

In this section, we present the qubit realization and 
give the effective Hamiltonian for a single hole in a quan- 
tum dot molecule (QDM) composed of vertically-stacked 
self-assembled InAs/GaAs quantum dots (QD) separated 
by a small tunnel layer. The detailed geometry of the 
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QDM, which we use for the present study, is shown in 
Fig. [lja). Two pyramidal dots (height h, width w) are 
stacked on top of each other and separated by a barrier 
of distance d. The QDM is exposed to a static magnetic 
field, as well as a time-dependent electric field applied 
along the [001] growth direction. It is controlled by a 
gate bias and used to modulate the hole g factor by shap- 
ing of the hole wave functionPSl Next to the interaction 
of the hole with the externally applied control fields, all 
contained in the effective Hamiltonian Hq, there will be 
additional (unwanted) interactions with the solid-state 
environment. Therefore, we account for the interaction 
of the host-lattice nuclear spins with the hole, develop a 
phenomenological description of the hole-phonon interac- 
tion, and give the resulting equation of motion on which 
our analysis is based. 

The first task is to compute the effective interaction 
of the hole with a constant external magnetic field. It 
is characterized by a g tensor, which depends on the ex- 
ternally applied electric field, as shown in Fig. |ljb) . 20 
We have performed three-dimensional 8-band envelope 
function calculations, including external fields, strain, 
and piezoelectric polarization, in order to determine the 
QDM heavy hole and light hole components of the ground 
|\&o) and first excited Zeeman state l^ i). De tails of this 
method have been published elsewhere! 21 * 22 ! Both states 
are predominantly heavy-hole (hh) like. However, light- 
hole (lh) contributions cannot be neglected. We write 
the hole wave functions as 



3,3* 

^{|^},i.e{44} 1 *e{0 >1 } 1 (1) 

where F( k 'd>i*) (E,r) denotes the envelope function as- 
sociated with the basis function ip^^^ir), which trans- 
forms like the eigenfunction \ j,j z ) of the angular momen- 
tum operator J, and f2 is the volume of the unit cell of 
the crystal. For zero electric field E and for vertical mag- 
netic field B — (0, 0, 10) mT, the hh and lh contributions, 
respectively, are given by 



dr\F {k ^' 2 ^l 2 \E,r)\ 2 » 0.915, 

J d?\F<- k & 2 >±V 2 \E,?)\ 2 « 0.074. 

We have chosen our coordinate system along the cu- 
bic axes, where z is the growth direction of the QDs. 
The conduction and split-off (SO) band contributions 
are neglected, since they are smaller than 1%. In 
Fig. [21(a) and[2lb), respectively, we present contour plots 
of \FV;V^-m( E ,r)\ 2 and \ F W>'-W{E, f)| 2 in the 
x = plane. The corresponding coordinate system is 
depicted in Fig. [T](a). 

The energy difference between |\E' 1 ) and the next higher 
energy eigenstate is larger than 1 meV, versus a 
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FIG. 1: (a) Vertically stacked pyramidal quantum dots, dot 
spacing d = 1.5 nm, height h — 2.5 nm, width w = 15 nm. 
The coordinates x, y and z denote the [100], [010] and [001] 
directions, respectively. The origin of the coordinate system 
is indicated by 0. The external electric field points along the 
[001] direction, i.e. E = (0,0, E). (b) Electric field depen- 
dence of the g tensor elements. The solid and the dashed line 
denote g 001 (E) and g 110 (E), respectively. (Data taken from 
Ref.EOl) 



splitting of ss 0.4 ^eV between the lowest two states. 
Hence, the system near the ground state is well de- 
scribed by a two-level system with basis states \^o) and 
l^i). These are essentially linear combinations of upper- 
and lower-dot hole states, with the admixture depending 
upon the value of the external electric field. The effective 
Hamiltonian of the pseudo-spin system readspHl 



H a = f -fa-g(E)-B, 



(2) 



where a, g{E) and -B, respectively, denote the Pauli 
matrix vector of the pseudo spin-^ system, the electri- 
cally tunable hole g tensor and the externally applied 
magnetic field. Note that Eq. |2]) is given in the basis 
, l^o)}- By choosing a constant magnetic field of 
the form B = (B, -B, B), Eq. (2| takes the form 
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Throughout this paper, we choose a value of B — 10 mT 
and an external electric field pointing along the [001] di- 
rection, E = (0,0, E) 



A. Hole Nuclear-Spin Interaction 

In addition to the externally applied magnetic field, 
the hole experiences an effective magnetic field that re- 
sults from the nuclear spins of the host lattice. This is a 
consequence of the non- vanishing dipole-dipole hyperfine 
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lh envelope function jF (1;3/2 '- 1/2) (£, f)| : 





FIG. 2: Spatial dependence of the envelope functions 
|F (1;3/2y '' ) (£;,r)| a within the x = nm plane [see Fig. Qa)] 
for (a) the hh contribution with j z = —3/2 and (b) for the 
lh contribution with j z = —1/2. All envelope functions are 
given for E = kVcm" 1 and B = (0,0, 10) mT. The cross 
sections of the pyramidal quantum dots are indicated by red 
dashed lines. 



interaction (in contrast to the Fermi contact hyperfine 
interaction that vanishes for wave functions of p symme- 
try), as well as the coupling of the h ole orbital angular 
momentum to the nuclear spinsP^3 We now determine 
this effective field. 

The corresponding interaction Hamiltonian of a single 
nuclear spin with a hole has the formp 3 -' 



Hi = 2n Bli P 



Pi 



Pi 



(4) 



where P denotes the i-th nuclear spin operator and 7, 
and hb, respectively, denote the gyromagnetic ratio of 



the i-th nuclear spin and the Bohr magneton. The hole 
spin operator is denoted by s, pi is the distance vector 
between the hole and the z-th nuclear spin (located at 

p — 0), and l % = pi x p denotes the hole orbital angular 
momentum operator. The matrix elements of Hi are now 
calculated in the {l^o) , l^i)} basis defined above. We 
closely follow the approach given in Refs. IT21 and IT41 and 
obtain 



3z>3* 



(5) 

where Ri denotes the position of the nuclear spin in 
the coordinate system of the envelope functions. The 
matrix elements with respect to the 4 basis functions 
j^(§ >D ; ■ 2) , <0(f>~§) ; ^(f'~f)| are given b}®' 
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where I± = P x ± P and P z are the nuclear spin operators 
and 



8p,Bjiti£l 



dp 



<P)\ 2 



(7) 



The integration in Eq. ([7| extends over the dominant 
part of the interaction defined by a radius Rq around the 
nuclear spin under consideration and n(p) is the radial 
part of the basis functions. The details of the calculation 
are given in Appendix [5] Finally, the interaction of the 
ensemble of nuclear spins with the two-level system can 
be cast into the form 



Hnuc(B n ) — H} = 



(8) 



where i runs over all N w 10 4 — 10 5 nuclear spins inter- 
acting with the hole. This has the form of the interac- 
tion of the pseudo spin with an effective operator- valued 
magnetic field B n . The dynamics of the hole is much 
faster than that of the nuclear spins. This allows one 
to employ a quasi-static approximation for B n for the 
time period of a single measurement (initialization, ma- 
nipulation and readout) of the hole spin qubit so th at B n 
may be approximated by a classical constant vector HH1241 
During the time it takes to perform 10 3 — 10 4 repe- 
titions of the measurement, the effective nuclear mag- 
netic field varies significantly, leading to inhomogeneous- 
broadening-type dynamics™ The simplest way to take 
into account this variation is to treat the vector com- 
ponents of B n — (B*, B^, B^) as random variables with 
Gaussian probability distributions) 3 -' 4 ^ 



P{B a )=P{Bl)P{Bl)P{Bl) 1 
with P(Bl) = cxp [-{B^ 2 /(2Ai 



(9) 
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Here, P{B l n ) is the probability of finding a value B l n for 
the effective nuclear magnetic field along the i-direction 
and Aj = (B^B^) — (B^) denotes the corresponding 
variance of the effective magnetic field fluctuation. Wc 
assume an "infinite temperature" nuclear-spin density 
matrix, i.e. p n = (21 + 1) N I. 4 Hence, the mean values 
vanish for all directions i. Furthermore, we assume 
that the spin bath is uncorrelated, i.e. = 0. 

Finally, the total Hamiltonian of the system reads 



H(B n ) — Hq + H nuc (B n ). 



(10) 



We note that H(B n ) is time-dependent via the electric 
field that is applied to control the hole dynamics. 

In a recent experiment, the dephasing time of a hole 
spin in a single quantum dot has been measured as 



-^2 h ~ 100 nsPSl This corresponds to a variance of A z 
h/(p, B T2 h ) ~ 0.1 mTpS We use this value for our simu- 
lation of the hole nuclear-spin interaction. The variances 
A x and A y are calculated in Appendix [Bj We find that 
A x /A z « A y /A z « 10" 1 . 



B. Hole Phonon Interaction 

The coupling of the hole to acoustic phonons via the 
piezoelectric and deformation potential interaction leads 
to additional dephasing and relaxation!^! To account for 
these mechanisms, we employ a Lindblad modelP^ The 
dissipator reads 



+ -y- WzPv z - p] , 



with Lj,, and T* h denoting the relaxation rates for the 
transitions -> |* ), |*o) -> and the pure de- 
phasing rate, respectively. The braces denote anticom- 
mutators. For simplicity we set r* h = = T± = T. 
The range of reported relaxation times varies significantly 
with tem perature and the externally applied magnetic 
fieldPZES p or this work, we choose a relatively conser- 
vative value of T hh = 1 /is = l/(r t + rj = 1/T. We 
note that the hole relaxation time for quantum dots due 
to phonon int eraction in creases with decreasing external 
magnetic field! 11 ' " 5 ™ ! 



C. Hole Spin Dynamics 

Inspection of Fig. |T|(b) reveals that the g tensor com- 
ponent g 001 (E) cannot be tuned to 0, in contrast to 
g 110 (E). For a simpler description of the hole-spin dy- 
namics, it is therefore useful to switch to a rotating frame 
\ip) — > = U2 \ip) which rotates around the z-axis with 
a frequency given by 

=n B g mi {E*)B/h. (11) 



The transformation is characterized by the time depen- 
dent unitary operator J7 2 = exp [iu>*to~ z /2], where the 
electric field E* is defined by the relation 



g iw {E*) =0. 



(12) 



In addition, we perform another time independent rota- 
tion U\ = e _I s <Ti , which corresponds to the pseudo-spin 
rotation cr x — a y —¥ \[2a x in Eq. pj). In this rotating coor- 
dinate system (labeled by a tildej^ the Lindblad equation 
for the density matrix reads 



P{B n ,t) 
dp{B n ,t) 
dt 

HABt.) 



U 2 U lP (B u ,t)ulul 
-~ \H r (B n ),p(B n ,t) 
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U 2 U 1 H[B ri )Ulul 



dt 2 



(13) 
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p(B n ,t) 



V 



p(B n ,t) 



with H(B n ) given in Eq. (10 1. It can be shown that 



the form of the dissi- 



pator is invariant under the coordinate transformations 
described above. The density matrix has to be aver- 
aged over the effective nuclear magnetic field B n of each 
measurement. It is calculated by averaging over typi- 
cally M — 3000 values of the random effective nuclear 
magnetic field with probability distributions as given in 
Eq. (9), using 



~p{t) = {l/M)Y d P{Bn) 



(14) 



M 



III. RESULTS 
A. Pulse Shape Optimization 

In the previous section we have detailed the hole-spin 
qubit and identified a complete control mechanism. We 
are now in a position to find optimal electric fields that 
perform any type of qubit transformation. However, the 
dependence of the g tensor elements on the electric field E 
is a complex one, as can be seen in Fig.[TJb). Therefore, 
appropriate control fields can, in general, not be deter- 
mined analytically, particularly, if the finite rise time of 
the electric control is to be taken into account. Here, we 
apply optimal control theory in order to determine both 
realistic and simple pulse shapes. 

We start by outlining how to characterize qubit oper- 
ations. For single qubit systems, quantum gate transfor- 
mations can be described as rigid rotations of the Bloch 
sphere. For unitary dynamics, this rotation can be fully 
described by the time evolution propagator U(t) corre- 



sponding to the Hamiltonian H r (B n ) given in Eq (13 1. 
However, it is more advantageous to employ the so-called 
process tomography matrix (PTM) x(t). For strictly uni- 
tary dynamics, this matrix takes the simple formp^l 



X (t) = U(t)*®U(t), 



(15) 
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with U it)* denoting the complex conjugate of the matrix 
representation of U(t). The PTM eliminates physically 
irrelevant global phases of the propagator U — > e z *£7, 
can be readily obtained in experiments^* ^, and can be 
employed for non-unitary dynamics as well 1241251 

The next step is to define a cost functional that quan- 
titatively reflects how accurately a control field achieves 
a given unitary gate transformation. We use the PTM 
of Eq. ( 15 ) that implicitly depends on the electric field, 



x(t) = x\E(t); t], and seek electric pulse shapes that min- 
imize the following cost functional,^ 

J[E] = tr { [ X (tf) - X d] [X(tf) ~ Xd? } , (16) 

where xd denotes the ideal PTM and the interval (0, tf) 
is the time-span allowed for the gate transformation (gate 
operation time). A perfect implementation of the given 
unitary gate transformation corresponds to J = 0. Other 
distance measures for open quantum systems are given in 
Refs. l37ti4Tl We choose the following analytic form of the 
electric control field, 




t - £ At, 



1 + tanh 



(17) 



where Aj and Atj are the parameters to be optimized and 
E* is the working point defined by Eq. (12). The pulse 



form of Eq. (17 1 corresponds to a sequence of M voltage 



steps of amplitude Aj , each of duration Atj . The finite 
rise time of the pulses is determined by the parameter a 
that we set to a = 8.79 ns _1 . The optimal control field 
E opt (t) is obtained from the minimization procedure 

Jopt[£o P t] = r min J[E(Aj, Atj)}. (18) 

Thus, a value J > reflects suboptimal pulse shaping 
and/or the presence of decoherence and relaxation ef- 
fects. In addition to J, it is customary to define a fidelity 
loss of the gate operation by 

AF[E} = {J[E]/J max )i (19) 

with J max = 2n 2 , where n denotes the number of basis 
states of the quantum system, i.e. n = 2 for single qubits. 
A perfect execution of the gate operation corresponds to 
AF[E] = 0. For brevity, we will write AF instead of 
AF[E] throughout this work. 



We perform the cost functional minimization Eq. ( 18 ) 
for both the Hadamard operation (which transforms the 
axes of the Bloch sphere as follows: x — > z, z — ► x and 
y — > —y) and a 7r/2 rotation around the y-axis. The 
corresponding unitary evolution operators Uh and {/a., 
and the ideal PTM xh and x~i respectively, read [see 
Eq. pj] 



U H = 



XH 



1 

71 

(Uh) 



1 

-1 

U H 



Uk = 



1 



-i 
i 

•U: 



We set the gate operation time tf — 10 ns. The min- 
imization is executed within a differential evolution al- 
gorithmP2l For the Hadamard gate, the algorithm con- 
verges to an optimal control field that consists of three 
pulses, depicted in Fig. [3]ja). The n/2 rotation, on the 
other hand, is realized by an optimal control field con- 
sisting of four pulses, as shown in Fig. [3jc). Trajecto- 
ries of the Bloch vector for both transformations with 
equal initial states (^i) are given in Figs.gb) andgd). 
If we neglect decoherence and relaxation, the fidelity 
losses due to imperfect pulse shaping are of the order 
AF s» 0.001%. When decoherence and relaxation are 
included, the fidelity loss increases to AF s» 1%, which 
gives excellent performance. For the Hadamard gate, the 
real parts of the matrix elements (to, n) of the optimal 
PTM, Ke[x(tf) m n], as well as the deviation of the optimal 
PTM from the ideal PTM xh, are shown in Figs. Qa) 
and Kb), respectively. 
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FIG. 3: (a) Optimal electric pulse shape E opt (t) for the 
Hadamard transformation and (c) for the n/2 rotation. The 
dot-dashed red lines correspond to E = E* given in Eq. ( 12 1. 



(b) The trajectory of the Bloch vector for the Hadamard gate 
and (d) for the 7r/2 pulse. The final states of both transforma- 
tions coincide for the initial state $i, though the trajectories 
of the Bloch vectors differ. 



B. Hole Spin Echo 

It is frequently important to analyze the characteris- 
tics of dephasing processes of a quantum system that are 
not associated with inhomogeneous dephasing, such as 
caused by the hole nuclear-spin interaction discussed in 
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(a) 



(b) 



ReDKWmn] 



Abs[*(tr)-* H ] mn /1(r 





FIG. 4: (a) Real part of the matrix elements (m, n) of the 
optimal Hadamard PTM xiff) 0°) Deviation of the optimal 
PTM with respect to \h (relaxation and dephasing included). 
The corresponding fidelity loss is AF « 1%. 



Sec. (A) in the present case. This can be achieved with 
spin-echo experiments.^ The decay of the peak value 
of the spin echo recovery signal gives information about 
the additional cumulative dephasing rate and the time- 
dependence of the coherence loss (e.g. exponential versus 
polynomial decay). 

In Fig. [5^ a) , we show the crucial steps that are needed 
to perform a spin echo experiment for the hole-spin qubit. 

(i) At time t = 0, one applies a ir/2 pulse such 
as given in Fig. pic). It transforms the ini- 
tial qubit state |Wi) into the superposition state 

iA/2(|* ) + l^i))!^ 

(ii) Subsequently, the Bloch vector of the qubit evolves 
according to the system dynamics, including dissi- 
pation, decoherence and the effective magnetic field 
from the nuclear spins. The latter causes a rotation 
of the Bloch vector around a random axis, given by 
the direction of B D . 

(iii) At time t\ > one applies a tt pulse that is com- 
posed of two subsequent tt/2 pulses. It rotates the 
Bloch vector by an angle tt around the y-axis of the 
Bloch sphere. 

(iv) Next, the Bloch vector evolves again and 

(v) at time t = 2t\ the coherence, represented by 
Re[poi(i)], is partly restored. 

Here, Po\{t) is the off-diagonal matrix element of 
the density matrix p(t) in the pseudo-spin basis [see 
Eq. ( 14 1] and serves as a measure for coherence. For 
t\ ss 0.4 [is, a plot of its time evolution during a spin- 



echo experiment, as obtained within our model, is shown 
in Fig. [5jb) . Spin-echo signals corresponding to dif- 
ferent pulse separation times t\ are given in Fig. ^c). 
Since relaxation and pure dephasing due to the spin- 
phonon interaction are included in our simulation, the 
restoration of coherence is obtained as imperfect. The 
peak values of the echoes decrease approximately as 
~ exp |— 2ti (r t + T l )/2 + r* h | , indicated by the red 



dot-dashed line in Fig. [5jc). We note that the non- 
vanishing transversal nuclear magnetic field components 
B* and B\ can also lead to a reduction of the peak value 
of the echo signal. However, we analyzed this effect and, 
for the present system, found it to be negligible com- 
pared to both the longitudinal contributions B„ and the 
phonon interaction. 
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FIG. 5: (a) Schematic spin-echo sequence as described in the 
main text. The red (dark gray) arrow denotes the Bloch vec- 
tor, (b) Time dependence of the coherence Re[/5oi(t)]. The 
spike at t as 0.4 ^is is a consequence of the 7r-pulse. At 
t ~ 0.8 [is a partial revival of coherence can be observed, 
(c) The solid black lines denote spin-echo signals for differ- 
ent pulse separation times ti. The decrease of the peak echo 
signal is denoted by the red dot-dashed line. 



C. Pulse Timing Imperfections 

As mentioned above, the 7r-pulse of the proposed spin- 
echo sequence is composed of two subsequent tt/2 pulses. 
Due to the form of the Hamiltonian in Eq. (13), these 



composite pulses should be applied at special times (in- 
teger multiples of = 2tt/lo* after state initialization) 
in order to preserve the high fidelity of the gate trans- 
formation. This requirement is due to the effective mag- 
netic field g(E) ■ B along the z-direction that cannot be 
tuned to zero by electric means. The dependence of the 
fidelity loss AF with respect to a pulse delay time error 
S between two subsequent 7r/2-pulses is shown in Fig. [6] 
This figure illustrates that the additional fidelity loss AF 
is smaller than 1% as long as 5 < 10 ps, which should be 
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readily within reach of present day experiments. The sen- 
sitivity on the delay time error can be further decreased 
by reducing the magnitude of the external magnetic field. 

(a) (b) 




(n-l)tR n b> (n+Dtn 6[ns] 



FIG. 6: (a) The schematic ideal pulse sequence is denoted by 
the black solid line. The second imperfect it/2 pulse (red/gray 
dashed line) is delayed by S. (b) Double logarithmic plot of 
the fidelity loss AF vs. the delay time error 5. 



IV. SUMMARY 

We propose a qubit realization in form of the spin of 
a hole in a quantum dot molecule, which is controlled 
all-electrically by g tensor modulation, and evaluate its 
performance regarding controllability and dissipative ef- 
fects. An effective qubit model is derived from a detailed 
electronic structure calculation, as well as the inclusion 
of the interaction with nuclear spins of the host lattice 
and phonons. On its basis, we predict that high-fidelity 



gate operations for a single hole within a quantum dot 
molecule are experimentally feasible and promising. The 
qubit can be fully controlled by means of electric g tensor 
modulation. We use optimal control methods in order to 
determine simple electric control pulses, as illustrated for 
the Hadamard gate and for a ir/2 qubit rotation around 
the y-axis of the Bloch sphere. The performance of these 
gate transformations is tested with respect to dephas- 
ing and relaxation due to the interaction of the hole spin 
with the surrounding nuclear spins of the host lattice and 
due to the hole-phonon interaction. For electric pulses of 
a duration of 10 ns, we find that qubit manipulations 
can be performed with a remarkably low fidelity loss of 
AF ps 1%. In addition, we propose a spin-echo experi- 
ment that allows one to completely cancel the inhomo- 
geneous dephasing due to hole nuclear-spin interaction, 
as well as a study of additional dephasing mechanisms. 
We also investigate the influence of pulse-timing imper- 
fections on the gate fidelity. We find that the additional 
fidelity loss is < 1% for pulse delay time errors of less 
than 10 ps. 
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Appendix A: Matrix Elements of the Hole Nuclear-Spin Hamiltonian 



In this section, we derive the matrix elements of Eq. |5]). For simplicity, we omit the nuclear spin index i. The 
interaction Hamiltonian of a nuclear spin I (located at the origin) with a hole (located at p) readsj^ 



Hi(p,p) = 2pBjI 



I s 3p (s • p) 



P &mn ^PmPr, 

~JpF~ 



(Al) 



where we used the Einstein summation convention for the indices m and n. Here 7, s and I = pxp, respectively, denote 
the gyromagnetic ratio of the nuclear spin /, the hole spin and the hole orbital angular momentum. In Eq. (All, 
we explicitly denoted the spatial and momentum dependence of Hi. The hole wave function 'F can be written as a 
product of envelope functions F and angular-momentum basis functions ipvJ*) t 

^ k {^ = ^Y j F^\E,r)^\r), 3, € {±f»±^} > J € {§, ^} , k € {0, 1} , (A2) 



or, using spin-resolved zone center valence-band Bloch functions tp^°^ , as 

^ k (r) = VQj2 F{k;i ' a) ( E ^ (i ' a) ^ ie{X,Y,Z}, ae{t,|}, k e {0, 1} . 



(A3) 



The matrix that describes the transformation between t/>CjJ*) and ^(*- a ) is given in Eq. ( A7 ) 122143 The wave functions 
^o(r) and ^i(r) denote the ground and the first excited Zeeman state of the quantum dot molecule. We neglect 



states with higher energy (see Sec. [TTJ) and calculate the matrix elements of the interaction Hamiltonian, 

Hi(p,p) = 



(*„IW,p)|*i) (*o|#i(p,p)I*o) 



(A4) 



with 



where -R denotes the position of the nuclear spin with respect to the hole envelope function coordinate system. Since 
the hole nuclear-spin interaction is short ranged, we can assume that it is non vanishing only within a sphere of volume 
Vr and radius Rq around the nuclear spinP^ Furthermore, the envelope functions are approximately constant within 
this sphere, i.e., 

(%|#i(p,p)|^) « n [F (k *' a) (E,R)\* F^\E,R) J dp[^' a \p)\* Hiifr-ihV^'Viff), 



Vr 



£ [F^ k ^ (E, R)] * F^>V ( E , R)V ljaP , 



(A5) 



where we used the spatial periodicity of il)( l ' a '(p) and (r| H\(p,p) \p) ~ H\(p, —ifiS/p)5{p — f). In fact, it can be shown 
that Vr is well approximated by the unit cell around th e nuc lear spin under consideration. Contributions due to long 
range interactions lead to corrections of the order of 1 % IMiil Following Ref. [T^l we employ a spherical approximation 
of the basis functions il)( l,a \p). Using spherical coordinates (p, the corresponding approximate basis functions 
read, 



^ x > a \p,e^) 
^ (1>) (p,M) 
^ q) (p,M) 



k(p) sin (9) cos ((f)) \a) , 
n(p) sin (8) sin (4>) \a) , 
n(p) cos (9) \a) . 



Here, n(p) denotes the radial part of ip( l ' a \p,9,(f)) and |a) is a ket in the spin basis {f, !}• An estimation of the 
magnitude of the hole nuclear-spin interaction is given in Ref. 1131 where the authors used a linear combination of 
atomic orbitals for the hh Bloch function. In the basis ij;^-^ [a, (3 € {t ii} an d hj G {-^> ^ ancl by employing 
the aforementioned spherical approximation, the integral Vij a p [Eq. (A5)] can be expressed as, 



dp 



(A6) 



A discussion of modifications of the interaction at small p due to the finite size of the nucleus can be found in Ref. l45l 
Since we want to express the hole nuclear-spin interaction in terms of lh and hh contributions, we chan ge the basis 
representation of Vij a p from Vifj^ l ' a H to (ip^'i*)}. The transformation is governed by the unitary matrix E2144J 



1 



-v^ iV3 

0-10 i 2 

10 i 2 

iV3 

~V2 iV2 -V2 



(A7) 



-V2 -n/2 V2 
where the new matrix representation V is given by = Up, a Vc,\U\ u . Here, we use the shorthand notation 

P,f = (j,jz) G 



3 3 
2' + 2 



3 1 

2' + 2 



3 1 

2' _ 2 



3 3 
2'^2 



1 1 

2' + 2 



1 1 

2'^2 



a, A = (i, a) G {(X, t), (X, |) , (y t) , (y |), (Z, t) , (Z,i)} ■ 



(A8) 
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The transformed matrix V^ v in the 



basis thus reads, 



- 8^ B jhn r \k{p)\* 

»V = ~. / dp 





#- 






























-/« _ 





(A9) 



Split-off contributions (j — |) have been neglected as discussed in Sec.[IT] We finally arrive at the expression, 



for the {l^o) i l^i)} matrix representation of the hole nuclear-spin interaction [for the definition of Vj> j x see Eq. ( A9 )]. 

The Wigner-Eckart theorem can be used to determine the matrix elements of the hole nuclear-spin interaction as 
weilP^HUl In addition to the calculation presented above, we decomposed the spatial part of Eq. (Al) into spherical 
tensors L 



and Q t - 



l m /p 5 



It. 



P $mn ^PmPn 



-> Qr. 



which correspond to angular momenta j = 1 and j — 2, respectively. The Wigner-Eckart theorem can then be 
readily applied to obtain the matrix elements in terms of Clebsch-Gordan coefficients and so-called reduced matrix 
elements. However, in order to get a relation between the reduced matrix elements of the spherical tensors l m and 
Qmn, approximations have to be employed as well. Hence, we decided to present the calculation of the matrix elements 
in terms of the spherical approximation given above. 



Appendix B: Variances of the Effective Nuclear Magnetic Field 



Here, we derive the expression for the effective nuclear magnetic field B n given in Eq. pi), as well as the ratio of 
the variances and A y with respect to A z [Eqs. pi]. We rewrite Eq. (AlOl into, 



(*k\H I {p,p)\V l ) = c A{E,R)I+ + A\E,Ii)I_+A z (E,R)I z 



(Bl) 



with A(E,R), A^(E,R) and A Z (E,R) denoting 2x2 matrices and with c given in Eq. Q. The corresponding matrix 
elements read, 



A k i(E,R) 
AUE,R) 



{V3 



1 

3 

\{ 

F^'i'^iR) 



F^i-i\R) 
F( fc -f<-i)(i?) 
p(i 



p(i; 



>(J2) 



pO 



3 

2 ■ 2 



l\R) - \F^'i'-i\R)] F^'i'-i\R)+ 



pi 1 ' 



where the electric field dependence of the envelope functions is omitted for brevity. 

For the given dot size, the hole spin interacts with typically N w 10 4 — 10 5 nuclear spins at positions Ri (i = 1...N). 
The full interaction Hamiltonian thus reads (with nuclear spin index i), 



H nuc = J2 H ^ = J2 Ci { [ A{ - E > + A ^ E > 



R/ 



A(E, Ri) — A\E, Ri) 



(E,R)li}. 



(B2) 



It can be cast into a pseudo spin-| form, where the hole interacts with an effective nuclear magnetic field B n , i.e., 
Hnuc — ' Bn- Contributions proportional to the identity operator of the pseudo spin- ~ space do not change the 
dynamics and have been omitted. The fc-th component of B n is calculated by projecting Eq. (B2| onto the Pauli 
matrix Ok (k £ {x, y, z}), 



(B n )k = — tr {a k H nuc } = — V \g^(E, R t )T l x + g*>(E, + g k n z {E, Ri)P z 

PB PB ; L 



(B3) 
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with 



g k n x {E,Ri) ee 2Rc [tr [a k A(E, Jfc)}] , g* y (E,Ri) = -21m [tr {^(B.ifc)}] , g*(E,Rt) = tr {a k A z (E, 



As noted in Sec. II A we employ a quasi-static approximation for the dynamics of the nuclear spins. We further 
assume that the spins are uncorrelated, e.g. (i^I™) — for j ^ m, and that only a single nuclear spin species with 
spin / is present. The "infinite-temperature" density matrix of the nuclear spin ensemble reads p nuc — 1(2/ + 1)~ N . 
Hence, the mean values (l x ), (l y ) and (/*) vanish and the variances of the nuclear spin components are (l x I x ) — 
(Iyly) = (PJl) = 1(1 + l)/3. Furthermore, we can set Cj = c. The variances thus read 



(B k n B*) 



Ms/ 



1(1+1) 



E 



[g k n x (E, 



gl v {E,Ri 



9n Z (E,Ri) 



(B4) 



The sum in Eq. (B4) is performed by a spatial sampling (10 4 points) of the lh and hh envelope functions. The value 
of the ratio of the variances with respect to A z is A x /A z sa A y /A z w 0.12. We then calculate A x and A y by 
setting A z = 0.1 mT (this value corresponds to the experimentally determined hole dephasing time T% h ss 100 ns, 
see Sec. II A). We note that the effective nuclear magnetic field and the variances depend on the externally applied 
electric field via the electric-field dependence of the envelope functions. The values for A k / A z given here have been 
calculated for E = kVcm" 1 and B = (0, 0, 10) mT. 
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